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Abstract 

We investigate the general asymptotic behaviour of Friedman- Robertson- 
Walker (FRW) models with an inflaton field, scalar-tensor FRW cosmo- 
logical models and diagonal Bianchi-IX models by means of Liapunov's 
method. This method provides information not only about the asymptotic 
stability of a given equilibrium point but also about its basin of attraction. 
This cannot be obtained by the usual methods found in the literature, such 
as linear stability analysis or first order perturbation techniques. Moreover, 
Liapunov's method is also applicable to non-autonomous systems. We use 
this advantadge to investigate the mechanism of reheating for the inflaton 
field in FRW models. 

1 Introduction 

A major difficulty in analyzing general relativistic problems, namely in cosmology 
and black hole physics, stems from the fact that the relevant field equations are 
nonlinear. This seriously limits the possibility of obtaining exact solutions, and 
makes it difficult to assess the degree of generality of the behaviour and special 
features of those exact solutions which are known. On the other hand, there 
is an increasing realization of the importance of the asymptotic behaviour of 
models, which provides the relevant features to be compared with the physical 
data available for each era of the universe. 

As a consequence, in recent years the emphasis in the study of cosmological 
models has drifted from the search for exact solutions to the analysis of the 
qualitative properties of the equations and of the long term behaviour of the 
solutions |1|]. The main tools used in this context have been perturbation methods 
and qualitative theory, especially linear stability analysis. Liapunov's method has 
been used in relatively few instances 0, |3| as a means to prove stability. 
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In this paper, we argue in favour of the use of Liapunov's method for the 
stability analysis of the possible attractors in a large class of cosmological models, 
as an alternative to linear approximation. One main advantage of the method is 
that it provides information not only about stability but also about the basins of 
attraction. Moreover, Liapunov's method is also applicable to non-autonomous 
systems j|, |fl. 

In its stronger version, Liapunov's method depends on the construction of 
a phase function F with the properties of being strictly decreasing along the 
orbits and having a minimum at the equilibrium point. Although stated in most 
introductory books in dynamical systems, this method is seldom used in the 
applications because it is difficult in general to find the form of the function 
F . However, there are some simple models in physics which indicate an obvious 
candidate for F such as in the case of the classical textbook application, the 
damped oscillator: the energy of the system is a positive definite form, and 
decreases along the orbits, thus verifying the condition of applicability of the 
method in its stronger version. 

A wide class of cosmological models share with the latter example the same 
form and qualitative dynamics. The equations of motion for all these models can 
be generically written as 

a 2 N 1. 

3 - = Enb* + M(a,b 1 ,---,b N ) (1) 
a i=0 A 

k = -3-bi- ^(o, b N ), i = l,2,...,N (2) 
a obi 

where a is non-negative. 

The degrees of freedom can be divided into two sets. The first corresponds to 
the global scale factor a which has, in the simpler models, a monotone behaviour 
in time, and the second set is formed by the quantities bi, i = 1, 2 . . . , N, which 
in many cases behave as damped oscillators and are thus especially well suited 
to the study and analysis through the construction of Liapunov functions. 

The paper is structured as follows. In Section 0, we briefly recall the basic 
results underlying Liapunov's method. In Section |3] we shall consider various 
examples for which a Liapunov function may be constructed that coincides es- 
sentially with the energy of the damped oscillator modes. For each case, the 
relevant dynamical information provided by the method is presented. First, 
in Section |3.1| , we study the basin of attraction of the de Sitter behaviour in 
Friedman-Robertson- Walker models with an inflaton field having a non-vanishing 
vacuum. We also assess the conditions for the existence of parametric resonance 



when the scalar field equation exhibits a periodic forcing term. In Section 3.2 



we consider the specific model of single-field reheating || [7|, [8], ^, |TI| and show 
that the results derived on the parametric resonance define the admissible class 
of potentials. Next, in Section |3.3|, we consider the convergence of scalar-tensor 
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gravity theories to general relativity in the case of Friedman-Robertson- Walker 
models |JTT], [12|, Subsequently, in Section |3T4] , we restrict ourselves to the 



simple case of diagonal Bianchi type-IX cosmologies with no matter and a posi- 
tive cosmological constant and study the validity of the cosmic no-hair theorem 
for this model. We shall recover in a simple manner, for this model, the result 
0. 



proved in 



Finally, the new information provided in each case by Liapunov's 



method is summed up in Section 



2 Reminder of the basic theorems 

Let us start by stating some general definitions. 

Given a smooth dynamical system x = f(x), x G IR n , and an equilibrium 
point xq, we say that a continuous function F : IR ra — > IR, in a neighborhood U of 
xq, is a Liapunov function for the point xq if 

(i) F is differentiable in U\{xo}, 

(ii) F(x) > F(x ) and 

(iii) F(x) < for every x G U\{x }. 

The existence of a Liapunov function F guarantees the stability of Xq] further- 
more, if the strict inequality F(x) < holds in U\{xq} then xq is asymptotically 
stable (see for instance 0) 

The neighborhood U is essentially constrained by condition (iii). So, when (iii) 
holds strictly, the method provides information not only about the asymptotically 
stability of the equilibrium point but also about its basin of attraction, which 
must contain the set U. This kind of information cannot be obtained neither by 
linear stability analysis nor by first order perturbation theory. 

Again in the case when (iii) holds strictly, the existence of a Liapunov function 
for the equilibrium point Xq has important consequences for the behaviour of 
time dependent perturbations of the above autonomous system. Consider the 
perturbed system 

x = f{x) +g(x,t), (3) 

where g(x,t) is smooth and bounded. Malkin's theorem [[| states that the so- 
lutions x(t; x(to),t ) of d||) will remain for all time inside a prescribed neighbor- 
hood B e of Xq, provided that we pick initial conditions x(to) sufficiently close to 
xq and that the magnitude of the perturbation ||<7(a;, t)\\ is sufficiently small for 
all t > to and all x G B e . This theorem guarantees the stability of xq for any 
sufficiently small time dependent perturbation, but it does not imply that the 
solution x(t; x(to),to) tends to xq when t — > oo, since it does not even require 
that g(x, t) should vanish as t — > oo. If this is the case, then a related result also 
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due to Malkin (see [[|) shows that an asymptotic property persists, in the sense 
that, as time increases, the neighborhood B e may be taken progressively smaller. 

Liapunov's method can also be used in the more general setting of time de- 
pendent systems which are not perturbations of autonomous systems (see [|J). 
Given a smooth non-autonomous system x = f(x,t), x G IR n , t G IR, and an 
equilibrium point x such that f(x ,t) = 0, if there exists a smooth function 
F : U C IR™ x IR — > IR, U a neighborhood of x , such that F(x , t) = F = const, 
F(-,t) — F is positive definite on U for every t G IR and F(x,t) < on U x IR 
then xq is stable. 

This result is analogous to the weak version of the theorem for the autonomous 
case. In the non-autonomous case, however, the conditions for asymptotic sta- 
bility are more restrictive. The point Xq is asymptotically stable if, moreover, 
F(x, t) is bounded from above, for all t, by a continuous increasing function of 
the distance \\x — xq\\, and F(x,t) is negative definite in U for all t G IR. 



3 Stability analysis 

3.1 The basin of attraction of de Sitter inflation 

During inflation the leading contribution to the energy momentum tensor is given 
by the inflaton scalar field <fi. The evolution of the Friedman-Robertson- Walker 
universes is described by 

where a is the (positive) scale factor of the universe, and k G { — 1,0, 1} distin- 
guishes the various spatial curvature cases (notice also that we adopt units in 
which 87rG = 1 = c throughout the paper). Under some assumptions on the 
scalar field potential V and on the initial conditions 0o> 0o an d cto; there exist 
solutions where the friction term 3H<p dominates in (|J) over and the potential 
term in (H) dominates over the kinetic term [[15], |I6 |. This is the inflationary 



stage, where the universe evolves with accelerated expansion. For definiteness, 
we shall consider the simplest model with V(</>) = A + m 2 2 /2, A > 0. 

With the change of variable H = a/a and b = 1/a equations @ and (|5|) 
become 

3H 2 = -3kb 2 + + V((f>) (6) 

b = -Hb (7) 
dV 

4 = -3H<j>-— (0). (8) 
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For future use we shall also need 

H = kb 2 - ^0 2 . (9) 

Our first goal in this section will be to apply Liapunov's method for au- 
tonomous systems to the characterization of the late time behaviour of the model 
associated to equations @, (|7|) and (§). This system has an equilibrium point 
at (0 = 0, = 0, b = 0) corresponding to a = oo where H = y\/3. Since 
H = a/ a, this equilibrium point corresponds to an infinite universe with expo- 
nential asymptotic expansion rate whenever A ^ 0. We shall therefore call the 
point (0 = 0, = 0, b = 0) the de Sitter spacetime. 

Let us prove that for k = —1,0 any initially expanding model will approach 
the de Sitter spacetime asymptotically for any O , O an d b , and that for k = 1 
the basin of attraction of the de Sitter attractor is given by b < JX/3. 

From (|G[), for an initially expanding model, H can be written as a function of 
b, and as 

H(4>, 0, b) = -L (skb 2 + ^ 2 + \ + ^m 2 2 ) V2 . (10) 

Notice that the argument of the square root is always positive for k non-positive. 
For th e cas e k = 1 this is also true if the initial value of the scale factor satisfies 
bo < \JX/3, because then b is a monotonically decreasing function of time. 

The asymptotic stability of the de Sitter attractor is obtained, as in the case 
of the classic damped oscillator, by using a Liapunov function given by 

F(6,0,0) = i 2 + i0 2 + K(0)>O. (11) 

From equations (0), (||) and (0) we have 

F(b, 0, 0) = -3H (0 2 + b 2 ) < 0, (12) 

since H((j), 0, b) > 0. Then, F given by (|IT|) satisfies (i), (ii) and (iii) for equations 
©j ®> an d (hi) holds strictly. This means that for k = —1,0 and A > the 
de Sitter spacetime is a global attractor, while for k = 1 and A > its basin of 
attraction contains all initial conditions (0 O , <fio, such that b < \JX/3. 
Consider now the following modification of the equations of motion 

3H 2 = -3kb 2 + ^0 2 + V{(j>) (13) 

b = -Hb (14) 
= -3#0-m 2 (l + £/(£)) 0, (15) 
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where / is a bounded, continuous, periodic function with period close to the free 
oscillating period of the unperturbated equation fll5l). 

We shall take equations (0), (PD , (|I5|) , as a toy model for the study of 
the problem of reheating, or parametric resonance, and defer to the following 
subsection its application to a specific model. 

Because (iii) holds strictly Malkin's theorem can be applied to this case. It 
states that the solutions of fll3"[ ), fll4]) and (|15|), will remain for all time inside 
a prescribed neighborhood of the infinity attractor provided that we pick initial 
conditions sufficiently close to it and that the magnitude of the perturbation is 
sufficiently small for all t. This theorem guarantees the stability of the infinity 
attractor for any sufficiently small time-dependent perturbation, but it does not 
imply that the solution tends to it as time goes to infinity. Nevertheless, we can 
look for some hypothesis concerning H which will entail the asymptotic stability 
of the equilibrium point = = 6 = 0. Let 

F{b, 0, 0) = \tf + \ (3a# (0, 0, b) + m 2 ) 2 + «00 + h 2 , (16) 
which is positive definite for < a < m. We obtain, using fllBP, (|T4D, (|I5|), 



F = - (3H - a) 2 - a (m 2 - + m 2 ef(t)^j 2 -m 2 £/(t)00 

- 3Hb 2 . (17) 
F will be negative definite if 

(3H - a) [m 2 - |# + m 2 ef(tj) - i m V/ 2 (t) > 0. (18) 

Since H{4> = 0, = 0, b = 0) = 0, this condition will hold for t large enough and 
e small enough whenever lim H(t) > a/3. 

t — >+oo 

There are two basic conclusions that can be obtained from these facts. First, 
parametric resonance can only occur when (|TJ) does not hold, that is, parametric 
resonance can only occur if H tends to zero as time goes to infinity, because, 
otherwise, the Liapunov function (|T6|), ensures that the equilibrium point = = 
&o = is asymptotically stable. This implies, for this toy model, that reheating 
can only be obtained from a potential that has a null minimum V(0) = A = 0. 
Second, even if H tends to zero as time goes to infinity, parametric resonance 
can only be observed if e is greater than some critical value because otherwise 
Malkin's theorem implies the stability of the infinity attractor. 

3.2 Single field reheating 

We now examine a specific model of reheating after inflation. According to the 
inflationary theory, almost all elementary particles populating the universe were 
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created during the process of reheating after inflation ]/]]. Reheating occurs due 
to particle production by the oscillating scalar field <fi. In the simplest models, 
this field (f) is the same inflaton field that drives inflation at the early stages of the 
evolution of the universe. After inflation, the inflaton oscillates near the minimum 
of its potential and this induces the quantum process of creation of particles |7J]. 
At the root of this process is a phenomenon of parametric resonance which gives 
rise to the excitation of certain modes of the quantum fluctuations of fields |17|] . 



We shall consider a simple model with the scalar field potential of the form 
V((f)) = Vo0", with Vo > 0, characteristic of the chaotic inflation scenario, and 
of the original studies of reheating. Our purpose is to ilustrate the possibilities 
underlying the aplication of the results of the previous subsection to this model. 
We aim at determining under what conditions does the mechanism of resonance 
occur, avoiding the drastic assumption often made in the analytic treatment found 
in the literature which consists in neglecting the Hubble rate of expansion || 
[7], |9|, [Hj. Furthermore, we also wish to show that the theorem on parametric 



resonance requiring the damping term to be vanishing in time and the amplitude 
of the periodic perturbation in the scalar field equation to be bounded from below 
permits to constrain the form of the potential in the case under consideration. 

We consider a FRW metric with non-positive curvature and shall focus on 
perturbations of the single scalar field ignoring the accompanying metric fluctu- 
ations. This is well justified on scales smaller than the Hubble distance, that 
is, satisfying p > aH, where p is the wave number of the perturbation, since 
relativistic effects play a lesser role on these scales On scales larger than 



the Hubble radius this approximation is more restrictive and, in particular, it has 
been recently shown to miss the resonant behaviour of the metric perturbations 
[|T^ , p0| , . With this caveat we shall thus consider the set of equations of motion 

3H 2 = -3kb 2 + ^<p 2 + V(<P) (19) 

b = -Hb (20) 
4> = -3H<f> - nVo^- 1 + b 2 Vl(p, (21) 

with n > 1 and x the spatial coordinates. 



Following fllOfl , we decompose the inflaton field into the sum of a classical 
inflaton field (p and a quantum fluctuation <50 in the form (f)(x, t) = (p(t) + 5(j)(x, t). 
We next concentrate in a specific mode of the fluctuation 5(f). This yields, keeping 



the first order terms in 5(f) p , 



(p + 3H (p + nVo^ 1 = (22) 
3H5<j) p + \n(n - 1) V^ n ~ 2 + p 2 b 2 } 5<j> p = 0. (23) 



Notice that n should be strictly greater than two for the existence of resonant 
behaviour. 
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Let dr/dt = b and denote the derivative d/dr by a prime. The equations of 
motion become 



3H 2 



b' 



V 3 



-3kb 2 + -b 2 V ' 2 + V(<p) 
-H 

2%ip' - nb^Vo^- 1 
b 

2 ^P 



n{n- l)b~ 2 V ip n - z +p 



(24) 
(25) 
(26) 

(27) 



with n > 2. 

Let T = (p 2 /2 = b 2 (p' 2 /2 and 7 = 2T/(T+V). Since 2T = nV holds asymptot- 
ically |22" |, we have 7 = 2n/(n + 2) and a simple argument shows that if ~ 6 3 ^ 2 
if A; = 0, H ~ b ii k = — 1. Consequently 6'/6 and b~ 2 (p n ~ 2 approach the infinity 
attractor as 



b'/b 



n-2u-2 



& (2n-2)/(n+2) jf ^ = Q 

6 if As = — 1 

6 (4n-16)/(n+2)_ 



(28) 
(29) 



From the previous discussion of our toy model we know that, in order to 



have parametric resonance, the damping term in (|27j) must tend to zero and the 
perturbing term in (^) must be bounded away from zero. The first condition is 
always satisfied, in view of (^) and of the fact that n must be greater than 2. 
From (^), the second condition is satisfied provided that n < 4. Notice that for 
n = 4 the term ip n ~ 2 b~ 2 in ( |27D oscillates with constant amplitude for all time. 
This is the case usually studied in the literature (see for example 0). However, 
( |2"9"| ) shows that reheating should be more effective for potentials with exponents 
closer to 2 in the allowed interval (2, 4]. 



3.3 General Relativity as a Cosmological attractor 
of Scalar- Tensor Gravity Theories 



General scalar-tensor gravity theories are based on the Lagrangian [11, O, 13 



L^ = ^R- -±-L g ^$ ^ v + 2C/($) + 16vrL m (30) 



where R is the Ricci curvature scalar of a space time endowed with a metric 
g^v, $ is a scalar field, u(Q>) is a dimensionless coupling function, U(Q>) is a 
function of and L m is the Lagrangian for the matter fields. They provide a 
most natural generalization of Einstein's general relativity, and their investigation 
enables a generic model-independent approach to the main features and cosmo- 
logical implications of the unifications schemes. The archetypal of these theories 
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is Brans-Dicke theory |23], 22] for which u;(<&) is a constant and the most distinc- 
tive feature of the theories defined by (|30|) is the variation of the gravitational 
constant, since G = 

Consider the homogeneous and isotropic Friedman- Robert son- Walker (FRW) 
universes and assume that the matter content of the universe is a perfect fluid 
with the usual equation of state p — (7 — l)p, where < 7 < 2. Here p and p 
are, respectively, the pressure and the energy-density measured by a comoving 
observer. Performing the conformal transformation of the metric El, E5| g a b = 



($/$*) g a b, where is a constant which normalizes the gravitational constant, 
to the so-called Einstein frame which corresponds to replacing the dynamical 
variables (a,t, $) with fyj^a, \f§~^i V 9 ) where 



rfln$_ / 16tt 
dtp VMMp) + 3)' 1 } 

the field equations reduce to the simple form 

3^ + 3^ = \<? + a-^M(<p) + V(<p) (32) 
a 1 a 1 2 

^ + 3-0 = -a-^M\v) - ^(ip), (33) 
a dtp 

where V(<p) = $*U ($(^)) / (8tt$ 2 (^)), M(^) = /i ($(^)/$,)^ 2+37/2 and p is a 
constant which fixes the initial conditions for the energy-density of the perfect 
fluid. We shall assume that M(tp) > (which corresponds to having G > 0) 
and that V(ip) = A + \m 2 ip 2 1 i.e., that V((p) has a non-degenerate positive global 
minimum at the origin. 

Let H = a/a and b = 1/a. It is easy to see that the former system is equivalent 

to 

3H 2 = -3kb 2 + -ip 2 + 6 37 M(v?) + V(ip) (34) 

b = -Hb (35) 

(p = -3Hip-b 3 m'(ip) - ^-(<p). (36) 

dip 

These equations define a dynamical system of dimension 3. We are interested 
in studying the stability of the equilibrium point at a = +00 corresponding to 



b = 0, H = yX/3 and <p — (p — 0, which is the de Sitter solution [pf] . 

Let us prove that for k = —1,0 any initially expanding model will approach the 
de Sitter spacetime asymptotically for any <p , ip Q and 60, and that for k — 1 the 



basin of attraction of the de Sitter attractor is given by b < y A/3. Consider the 
initial conditions H Q > associated with an expanding model, and b , (p and ip 
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arbitrary. Using the hypothesis M(tp) > 0, the expanding branch H = H(ip, ip, b) 
of (|35|) is given by 

H =-J§ {~ 3kb2 + \^ + h ^ M ^) + V (V)) 1/2 > • (37) 



As for the case k — 1, we take again Hq > and bo < yA/3 and, since, 
from fl35f ), 6 is a strictly decreasing function of time, it is possible for these initial 
conditions to write H = H{tp, 0, b) as in (pTf). 

The function defined by 

F(6, <p, <p) = ]^p 2 + A + V + 6 3 ^M(<^) (38) 

is a Liapunov function for the equilibrium point (6, 0) = (0, 0, 0). Indeed using 
(H), © and (|3), it is easy to see that 

F(b, <p, 0) = -3H (f + 7 6 37 M(^)) < , (39) 

and it is also straightforward to check that (i) and (ii) hold for F. Thus this 
proves that the equilibrium point (6, (p, <p) = (0, 0, 0) is asymptotically stable and 
shows that an attractor mechanism indeed exists when the matter content of the 
universe is a perfect fluid and the potential has a non-degenerate non-negative 
minimum. This means that for k = —1, and A > the de Sitter spacetime is a 
global attractor, while for k = 1 and A > its basin of attraction contains all ini- 



tial conditions ((p ,(p ,b ) such that b < yA/3. A similar attracting mechanism 
should also exist in the non-isotropic case. 

In general, linear stability analysis only provides information about the stabil- 
ity and asymptotic behaviour in a undetermined neighborhood of the equilibrium 
point. This can be easily seen in the following example. Consider the equations 
of motion d35|) , (|36D for 7 = 4/3 (for 7 ^ 4/3 the linear stability method can 
not be applied without additional hypothesis on M(ip)). The eigenvalues asso- 



ciated with the de Sitter attractor are A& = — J A/3, A^ = — V3A — \/3A — 4m 2 

and A^ = — v^A + V3A — 4m 2 . The linear approximation confirms the result 
obtained by the Liapunov method, the real part of the eigenvalues is always neg- 
ative, and shows that there exist two basic behaviours for the field associated to 
the monotone approach of the scale factor to the de Sitter attractor: an oscilla- 
tory mode (damped oscillations) if < A < Am 2 / 3 and a monotone behaviour 
(over-damped) for the field if A > 4m 2 /3 which corresponds, respectively, to a 
stable focus and a stable node in the plane (<p, if, b = 0). No information is ob- 
tained concerning the size of the basin of attraction for the de Sitter attractor 
with this method. 

Also notice that when 7 < 1/3 the second member of equation (36) is not 
differentiable at b = 0, which precludes the application of linear stability anal- 
ysis in its simplest form. This fact emphasizes the advantage of the use of the 
Liapunov's method. 
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3.4 Asymptotic behaviour for diagonal Bianchi type-IX 
cosmologies in the presence of a positive cosmological 
constant without matter. 

By definition, the general Bianchi type-IX spacetime has topology IR x S 3 , with 
a simply transitive action of the isometry group SU{2) on S 3 spatial slices. The 
metric of a general Bianchi type-IX model can be put in the form |26| 

3 

ds 2 = -dt 2 + a 2 it) Y \e p ] dw l dJ , (40) 

where uj % , i = 1, 2, 3, are isometry invariant one-forms on the three-sphere, a is a 
scalar, and (3 is a traceless 3x3 matrix. Both a and /3 are functions of the proper 
time t only. 

For a diagonal spacetime, let /3j, i = 1,2,3, denote the diagonal elements of 
the matrix (3. Only two of these quantities are independent since j3\ + fa + @3 = 
on account of the tracelessness of (3. We choose the independent variables to be 



h - (41) 

b2 = -^j= (Pi - fa) ■ (42) 

The vacuum Einstein equations with the cosmological constant take the form 
M 12711 



3^= 1 -(b 2 l +bl)+- 2 V(b u b 2 ) + A (43) 



1 ot r 

b l + 3% + -— = 0, i = l,2 (44) 
a a z obj 



with 



V(bt, b 2 ) = -e-^ cos h(V2b 2 ) + ^e- 4 VI b i (45) 

+ I e 2 v / f 61 [cosh(2 v / 26 2 ) - 1 
Let H = a/a and b = 1/a. It is easy to see that the former system is equivalent 

to 

3H 2 = ^(b 2 1 + bl)+A + b 2 V(b u b 2 ) = (46) 
b = -Hb (47) 
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dV 

h + 3Hbi + b 2 — = 0, i = l,2. (48) 
These equations define a dynamical system of dimension 5. This system has two 



finite equilibrium points b = 0, H = J A/3 and b = 2yA/3, H = both with 

b\ = b 2 = b\ = b 2 = 0. We are interested in the stability of the equilibrium point 
6 = associated with the de Sitter attractor. 

Notice that near the minimum (61,62) = (0)0) t ne potential ([45|) takes the 
form 

^(61,62) = ~+(bi+bi)+o 3 (b 1 ,b 2 ) 

= -1 + vihM), (49) 

where v(0, 0) = 0. 

We take H > and 6 < 2y / A/3. Since, from (p|) , 6 is a strictly monotone 

decreasing function with time the expanding branch of H = H(b, 61, 62, 61, 62) can 
be written for these initial conditions as 

1 ( 1 o n 1/2 

H = ^ {- (6? + bl) + 6^(6!, 6 2 ) + A - |6 2 } > . (50) 

The Liapunov function for the equilibrium point (6, 61, 6 2 , 61, 62) = (0, 0, 0, 0, 0) 
will now be 

Fib, h, 6 2 , 61, 62) = i (6? + 6 2 ) + 6^(6:, 6 2 ) + > . (51) 
Indeed using the equations of motion we have 



-3H 



hj + b 2 2 + h 2 v(bi,b 2 ) + b 2 



< (52) 



and thus it is a simple matter to check that all the conditions (i), (ii) and (iii) 
are satisfied. Therefore the equilibrium point (6,61,62,61,62) = (0,0,0,0,0) is 
asymptotically stable and the basin of attraction for de Sitter spacetime contains 
all initial conditions ^6 , 61 (0), 6 2 (0), 6i(0), 62(0)) such that 6 < 2^J A/3. In this 
case the linear stability analysis reveals only that the eigenvalues of the de Sitter 
attractor are all negative. 



4 Summary and conclusions 

In this paper we have applied the Liapunov's method to the study of the general 
asymptotic behaviour of cosmological models, namely: FRW cosmologies with 
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inflaton field; scalar-tensor gravity theories in FRW; diagonal Bianchi-IX vacuum 
cosmologies in the presence of a positive cosmo logical constant. For these models 
we recover through this method the well known stability result for the de Sitter 
solution and we show that the method also provides information about its basin 
of attraction. 

Furthermore, we have applied Liapunov's method to another problem of a 
completely different nature, namely the analysis of the conditions for the existence 
of the parametric resonance mechanism underlying the reheating process. We 
have concluded that in order to have parametric resonance, the damping term 
in the quantum perturbation equation must tend to zero and the forcing term 
must be bounded away from zero. This implies that for the basic model with 
potentials of the form V(tp) = Vo(p n , reheating may occur only when n e (2,4]. 
As far as we know, albeit the simplicity of the reheating model considered, this 
conclusion is new and fully illustrates the extension to which Liapunov's method 
may be utilized to draw rigorous results in cosmology. 
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